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Abstract The kinetics of color changes in keyaki 
(Zelkova serrata Makino) and sugi ( Cryptomeria japonica 
D. Don) wood during heat treatment were examined. The 
color of wood specimens treated at 90, 120, 150, and 
180 °C was measured by an imaging spectrophotometer 
and expressed using CIELAB color parameters. At any 
treatment temperature, values for L* and A E* ab decreased 
and increased in both wood species, respectively, with 
increased treatment time. Changes in a* and b* varied 
depending on wood species and treatment temperature. The 
color changes were successfully analyzed using the kinetic 
approach applying time-temperature superposition 
method. This approach elucidated and accurately predicted 
color changes during heat treatment. 
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Introduction 

Heat treatment of wood, particularly at temperatures below 
the thermal decomposition point of each wood component, 
has been widely studied for various purposes, including 
wood modification, wood drying, and accelerated aging [1, 
2]. Among the various properties of wood, color change 
during heat treatment has been frequently examined to 
develop chemical-free coloration methods [3] or to prevent 
staining and discoloration of wood during heating [4, 5] 
because color is esthetically important and affects the com¬ 
mercial value of wood. Many previous studies examined the 
effects of factors affected by heating processes such as 
temperature, duration, and moisture. However, individual 
phenomena in specific treatment conditions were described; 
no comprehensive study has been conducted to examine the 
color change behavior of wood during heat treatment. Color 
has also been studied as a nondestructive criterion for 
monitoring changes in other wood properties of wood. Pre¬ 
vious studies investigated treatment duration and tempera¬ 
ture within a limited range and for only a few wood species. 
Thus, the potential of color as a criterion remains to be 
determined. Some researchers found clear relationships 
between color and changes in other properties such as weight 
and mechanical properties during heat treatment [6-8], 
whereas others asserted that color is not suitable as a criterion 
for wood evaluation purposes, or that assessment techniques 
must be improved because results obtained in previous 
studies were unclear or occasionally contradictory [9-11]. 

In our previous studies [12, 13], color changes in hinoki 
wood (Chamaecyparis obtusa Endl.) and cellulose filter 
paper as a major component of wood were successfully 
analyzed using a kinetic approach, which allowed precise 
prediction of color changes at any temperature and time 
within the experimental ranges. Furthermore, the obtained 
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results of hinoki wood were expanded to elucidate color 
changes during natural aging [14], leading to the conclusion 
that a uniform reaction (thermal oxidation) causes color 
changes during both natural aging and heat treatment. 

In the present study, color changes in keyaki (Zelkova 
serrata Makino) and sugi (Cryptomeria japonica D. Don) 
were comprehensively investigated. These typical exam¬ 
ples of hardwood and softwood, respectively, are widely 
used in wooden buildings in Japan, such as hinoki wood. 
Specimens were heat-treated at various treatment durations 
and temperatures. A kinetic approach was used to analyze 
the temperature dependency of color changes in these two 
species. The results were then compared with those of 
hinoki wood and cellulose filter paper to elucidate changes 
in the color of wood during heat treatment. 

Materials and methods 

Samples 

A 100-year-old keyaki and a 180-year-old sugi tree were 
used. The trees were harvested in Japan and the wood was 
dried for more than 7 years under ambient room temperature 
conditions. A sample block of homogenous grain for each 
species was carefully selected from the area near the outer¬ 
most part of the heartwood. The size of specimens was 
60 mm (longitudinal direction) x 10 mm (radial direc¬ 
tion) x 2 mm (tangential direction). Air-dried density and 
annual ring width of each species are summarized in Table 1 . 
The specimens were dried in an air-circulating oven and a 
vacuum oven at 60 °C for 24 h each. Then, the specimens 
were heated in an air-circulating oven at 90, 120, 150, and 
180 °C for durations ranging from 10 min to approximately 
1.4 years (Table 2). Heating times were selected on the 
assumption that a 10 °C increase in reaction temperature 
halves the reaction rate, to obtain a similar degree of reaction 
intensity despite the differences in temperature [15]. 

Color measurement 

Despite careful selection of samples, the density varied 
more widely and the annual ring width were larger than 


Table 1 The average air-dried density and annual ring width with 
their standard deviations in parentheses 


Sample 

Air dried density 

Annual ring 

n 


(g/cm 3 ) 

width (mm) 


Keyaki wood 

0.66 (0.02) 

4.40 (1.35) 

255 

Sugi wood 

0.36 (0.02) 

1.50 (0.31) 

243 


n is the number of specimens used to calculate the average value 


Table 2 Treatment temperature and time for keyaki and sugi wood 



90 °C 

120 °C 

150 °C 

180 °C 

Treatment time (h) 

12 

1 

40 min 




3 




85.3 

10.7 

1.33 

10 min 


171 

21.3 

2.67 

20 min 


256 

32 

4 

30 min 


384 

48 

6 



524 

64 

8 

1 





1.5 


1024 

124 

16 

2 


3072 

384 

48 

6 


6144 

768 

96 

12 


9266 





12288 

1536 

192 

24 



3136 

384 

48 



5440 


84 



7680 

960 

120 


those of hinoki wood used in our previous study [12]. To 
reduce within-specimen variations in color caused by less 
homogeneity in tissue structure, specimen color was mea¬ 
sured with a spectroscopic imaging system. The system 
consisted of an imaging spectrophotometer (ImSpector 
V8E, JFE Techno-Research Corporation, Tokyo, Japan), a 
monochrome camera with a charge-coupled device, xenon 
light source, stepper table, and controlling computer. This 
system provided scanned digital images of all specimens. 
Each pixel in the images contains spectroscopic intensity 
data, which can be converted to various color parameters 
using spectrum-analyzing software (ImSpector Spectrum 
Analyser, JFE Techno-Research Corporation, Tokyo, 
Japan). An area of 50 mm (longitudinal direction) x 8 mm 
(radial direction) in each specimen was scanned at a 
scanning resolution of approximately 0.19 mm/pixel 
(134 dpi). The color value of each specimen was calculated 
as the average of all pixels in the measured area, which 
reduced the variations in color. 

The CIELAB color parameters (L*, a*, and /?*) were 
used to express color changes. The differences in param¬ 
eters AL*, Aa*, and A b* and the total color differences 
A E* ab were calculated using the following formulae [16, 
17]: 

AL* = L* - L* 

A a*=a*-o* 0 (1) 

A b* =b* -b* 0 

A E* ab = Val * 2 + Aa* 2 + A b* 2 , (2) 
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Table 3 The average initial color values Lq, 
standard deviations in parentheses 

a* 0 , and b* 0 with their 

Sample 

T* 

* 

a 0 

h* 

n 

Keyaki wood 

59.69 (1.78) 

11.79 (0.88) 

31.96 (0.98) 

155 

Sugi wood 

56.25 (2.43) 

11.63 (1.02) 

21.27 (1.03) 

90 


n is the number of specimens used to calculate the average value 


where L* is the lightness and a* and b* are the color 
coordinates under any treatment condition, respectively, 
and Lq, <Zq, and are the corresponding initial values 
obtained as the average of untreated specimens. The values 
of Lq, < 2 q, and b$ are summarized in Table 3. Four to five 
specimens were tested under each treatment condition and 
their average values and standard deviations calculated. 

Kinetic analysis 

A chemical reaction in a uniform phase can be described by 
the Arrhenius equation: 

fc = Aexp^-J^ (3) 

where k is the reaction rate constant, A is the frequency 
factor, E a is the apparent activation energy, R is the gas 
constant, and T is the absolute temperature at which the 
reaction occurs. In general, the temperature dependence of 
the reaction rate is expressed as E a . When changes in 
material properties were measured at different T, the values 
of k were obtained from their reaction rates. By plotting the 
logarithm of k versus the reciprocal of T, the values of E a 
can be obtained from the slope of the regression line for the 
plots, which is called the Arrhenius plot [18]. Interpolating 
or extrapolating the regression line of the Arrhenius plot 
allows us to determine the reaction rate at any given 
temperature, which is useful not only in understanding the 
reaction, but also in predicting changes in the property 
under any conditions. 

To determine the activation energy more accurately, a 
time-temperature superposition (TTSP) method has fre¬ 
quently been used for polymeric materials [19-22]. An 
Arrhenius approach applying the TTSP method is briefly 
described as follows. When changes in material properties 
at different treatment temperatures are plotted as functions 
of logarithmic treatment time and the plotted curves show a 
similar shape irrespective of treatment temperature, the 
plotted curves can be superposed to a single curve at a fixed 
reference temperature T ref along the logarithmic time axis. 
This shift distance used to superpose the curves is called the 
time-temperature shift factor a T and is given as follows: 


Clj 1 


tT 

fief 



where fi e f is the treatment time at T ref , and t T is the time 
required to obtain the same response at treatment 
temperature T. Combining Eqs. (3) and (4) thus gives: 


aj = exp 






where both T and T ref are absolute temperatures. Plotting 
the logarithm of a T versus the reciprocal of T is another 
way to calculate E a . Ding et al. [23] compared the TTSP 
method with the conventional method in Eq. (3) by ana¬ 
lyzing several sets of data for cellulose degradation from 
literature. They concluded that the TTSP method is more 
accurate than the conventional one. In our previous study, 
the color changes of hinoki wood [12] and cellulose filter 
paper [13] were successfully analyzed using a combination 
of kinetic analysis and TTSP. 

Though the curves plotted as a function of logarithmic 
treatment time are similar in shape, these curves cannot 
always be completely superposed, namely, a single 
smooth curve cannot be obtained only by a horizontal 
shift of data along the logarithmic time axis. In such 
cases, in addition to a horizontal shift, the data set may 
be shifted vertically along the logarithmic property axis 
to obtain a smooth superposed curve. This shift distance 
is called vertical shift factor b T . Vertical shift has been 
used to explain the influences of temperature on the 
viscoelastic behaviors of thermorheologically complex 
materials [24-26]. Though application of b T in the con¬ 
text of color changes during heat treatment has not been 
attempted, the TTSP method was modified in the present 
study using b T formulated as follows: 



where b T is the vertical shift factor at treatment temperature 
T, P T the property at T, and P re f the property at T ref . 

To calculate a T and b T , a temperature of 180 °C was 
chosen as T ref . Nonlinear regressions were performed to 
find the best models that fitted well to changes in color 
parameters at 180 °C, expressed as a function of t. 

All parameters of the model functions and kinetics 
were estimated using a nonlinear iterative curve-fitting 
method. The coefficient of determination (R ) and the 
root mean-square error (RMSE) were used as criteria 
for the fitting of the tested models to the experimental 
data: 



Em (A ~ F ) 2 

£f=i (Yi - Y t ) 2 
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Keyaki wood 




Sugi wood 



Treatment time (h) 


Fig. 1 Color changes of specimens treated at 180 °C in order of 
increasing treatment time 


RSME = 


Ef=i (Y t - Y t ) 


N 



S\ 

where Y t and Y t are the measured and modeled values, 
respectively, Y is the average of the measured data, and 
N is the total number of observations. The higher the value 
of R~ and the lower the value of RSME, the better did the 
model fit the data. The FlexPro 8 software package (Wei- 
sang GmbH, Germany) was used for all data analysis. 


Results and discussion 

Color changes during heat treatment 

Figure 1 shows the images of representative specimens. 
The specimens turned darker and browner with increased 
treatment time. The color changed uniformly from the 
outside to the middle of the specimens. Figures 2 and 3 
show changes in L*, < 2 *, b*, and A E* b of keyaki and sugi 


Fig. 2 Changes in L*, < 2 *, b*, 
and A E* ab of keyaki wood during 
heat treatment as a function of 
time and temperature (error bar 
standard deviation) 
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Fig. 3 Changes in L*, < 2 *, b*, 
and A E* b of sugi wood during 
heat treatment as a function of 
time and temperature (error bar 
standard deviation) 

* 

-«J 


■k 





Fig. 4 Changes in a* and b* of 
keyaki wood superposed using 
only horizontal shift factor a T 
(<error bar standard deviation) 
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wood, respectively, according to the treatment time and 
temperature. The values of L* decreased in both keyaki and 
sugi wood, indicating that they became darker. This result 
corresponded to the increase in A E* ab , which was largely 
dominated by decreasing L*. Decreases in L* and increases 
in A E* b during heat treatment have been well documented 
in various wood species [1-12]. In addition, changes in a* 


and b* differed depending on the wood species. In keyaki 
wood, the value of a * initially decreased, then increased to 
a peak and decreased again, whereas the value of b* 
decreased constantly with a shoulder. In sugi wood, the 
value of a* decreased slightly, whereas the value of b* 
initially increased and then decreased. This increase and 
subsequent decrease in a* and b* have been reported for 
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some other wood species such as beech, spruce, and pine 
[6, 7, 27], whereas a monotonous decrease in a* and b* has 
been reported for pine and eucalypt wood [3]. No signifi¬ 
cant change was previously reported in an analysis of sugi 
wood [28]. This inconsistency may be explained by the fact 
that treatment duration was insufficient in these reports. 

Comparing the results among treatment temperatures in 
each species, changes in L* and A E* ab showed similar 
behavior irrespective of treatment temperature, whereas 


Table 4 Coefficients of determination (R ) and root mean-square 
errors (RMSE) of regression curves using Eq. (9) for color changes of 
keyaki and sugi wood at 180 °C 



Keyaki wood 

Sugi wood 


R 2 

RMSE 

R 2 

RMSE 

L* 

0.9938 

0.294 

0.9962 

0.634 

<2* 

0.9957 

0.721 

0.9142 

0.307 

b* 

0.9970 

0.335 

0.9740 

0.542 

A F* 

L * r 'ab 

0.9964 

0.743 

0.9967 

0.614 


changes in a* and b* were temperature dependent. As 
shown in Figs. 2 and 3, the lower the treatment tempera¬ 
ture, the higher the peak values of a * for keyaki and b * for 
sugi, indicating that treatment at lower temperatures 
enhanced the increase in these chromatic parameters. 
Furthermore, the lower the treatment temperature, the 
smaller the decreases in b* of keyaki and a* of sugi wood, 
indicating that heat treatment at lower temperatures mod¬ 
erated the reduction in these chromatic parameters. No 
such temperature dependencies were observed in hinoki 
wood and cellulose filter paper [12, 13]. Thus, the results of 
our studies on heat treatment in several wood species at a 
wide range of durations and temperatures revealed that 
color changes in complex ways depending on the duration 
of heat treatment, temperature, and wood species. The 
results would contribute to the simple assessment of wood 
quality and understanding the empirical knowledge about 
quality changes. Almost all the main components of wood, 
i.e., cellulose, hemicelluloses, lignin, extractives, and ash 
contents are related to the color of wood and its changes 
[29]. Our previous study showed that cellulose is a major 
component that contributes to color changes of hinoki 


Fig. 5 Changes in color 
parameters of keyaki wood: 
best-fit model [Eq. (9)] and 
superposed data using only 
horizontal shift factor a T for L* 
and A E* ab and using both a T and 
vertical shift factor b T for a * 
and b* (error bar standard 
deviation) 
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wood, while other components also have effects on color 
changes [13]. Extractives of each wood species also make 
the characteristic color [30-32] mainly expressed by 
chromatic parameters, a* and b*. The different behavior of 
different color parameters would be caused by chemical 
changes of each component. However, CIELAB color 
parameters do not necessarily reflect a particular chemical 
change in wood. To refer to chemical changes, qualitative 
and quantitative analysis of different wavelength spectra 
and chemical components will be needed. 

Kinetic analysis 

In keyaki and sugi wood, changes in L* and A E* b were 
analyzed using the TTSP method and by applying only the 
horizontal shift factor a T . As shown in representative data for 
a* and b* in Fig. 4, single smooth superposed curves were 
not obtained for a * and b * by shifting only along the loga¬ 
rithmic time axis, though overlap of the inflection points of 
each curve was observed. This implies that application of the 
shift factor, a T , only was insufficient for accurate prediction 
of color changes in keyaki and sugi wood and that temper¬ 
ature dependencies should also be considered. Therefore, the 
vertical shift factor b T was applied in addition to a T to analyze 
the changes in a* and b*. 


Based on the values of R ~ and RMSE, the following was 
selected as the best-fit function: 

Cp + Cif + C 2 t 2 + C 3 f 3 

1 + c 4 t + c 5 t 2 + C 6 t 3 ’ ^ ’ 

where Q (i = 0, ..., 6) is a constant. Note that this function 
was applied for numerical analysis of color changes; it may 
have no meaning in terms of the reaction mechanism. 

Table 4 summarizes the values of R and RMSE for 
Eq. (9) fitted to the data set of color changes at 180 °C, 
indicating that Eq. (9) fitted well to the data set. Figures 5 
and 6 show the well-superposed data with the regression 
curves of keyaki and sugi wood, respectively, using only a T 
for L* and A E* lb and both a T and b T for a* and b*. Figure 7 
shows the Arrhenius plots for L* of keyaki wood as an 
example. The regression lines of the Arrhenius plots for all 
color parameters showed good linearity (R > 0.980), 
suggesting that color changes at 90-180 °C were caused by 
an apparently similar reaction. The apparent activation 
energies (E a ) calculated from each color parameter were 
higher than those of hinoki wood and cellulose filter paper 
[12, 13], particularly for L* and A E* clb (Table 5). This result 
indicated that color change reactions of keyaki and sugi 
wood apparently have higher energy barriers than those of 
hinoki wood and cellulose. Considering sugi and hinoki are 


Fig. 6 Changes in color 
parameters of sugi wood: best- 
fit model [Eq. (9)] and 
superposed data using only 
horizontal shift factor a T for L* 
and A E* b and using both a T and 
vertical shift factor b T for a * 
and b* (error bar standard 
deviation) 
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softwoods and keyaki is hardwood, the present results 
differ from those of previous studies, which suggested that 
the thermal degradation of hardwood proceeds faster than 
that of softwood because some components of hardwood 
are more reactive than those of softwood [33, 34]. The 
differences in E a values among color parameters may 
therefore be interpreted as differences in reactivity of the 
chemical components responsible for each color parameter. 



Fig. 7 Arrhenius plots for L* of keyaki wood, regression line 
(R = 0.997), and calculated apparent activation energy (E a ) 


Table 5 Comparison of apparent activation energies (E a ) calculated 
using the modified Arrhenius equation [Eq. (5)] among several wood 
species [12] and cellulose filter paper [13] 



E a (kJ/mol) 



Keyaki 

wood 

Sugi 

wood 

Hinoki 
wood [12] 

Cellulose filter 
paper [13] 

L* 

178 

180 

117 

125 

< 3 * 

132 a 

135 a 

95 

124 

b * 

133 a 

129 a 

114 

118 

A F* 

^ab 

188 

123 

113 

120 


a Values of E a were calculated by applying both a T and b T for <3* and 
b* of keyaki and sugi wood 


These results implied that reactivity may differ according 
to species and physical properties. However, further 
investigations including chemical analysis are needed to 
elucidate the E a values obtained in this study. 

The relationships between b T and treatment tempera¬ 
ture for a* and b* are shown in Fig. 8. As the reference 
temperature was 180 °C, b T at 180 °C was 1.0. The val¬ 
ues of b T represent the difference in the extent of decrease 
or increase in a* or between treatment at 180 °C and 
that at other temperatures. The values of b T decreased 
with increasing temperature toward 1.0 at 180 °C in both 
a* and £>*. In Figs. 2, 3, and 8, specimens were more 
reddish and yellowish at lower temperatures and b T 
expressed its degree. The relationship between b T and 
temperature was linear from 120 to 180 °C, but b T 
seemed to be saturated at temperatures less than 120 °C. 
This result indicates that treatment at lower temperatures 
keep colorimetric parameters larger and it may have an 
upper limit. Studies examining color changes at lower 
temperatures than those used in this experiment will allow 
further evaluation of b T . 

In the present study, the factors affecting temperature 
dependency expressed by E a and b T of color changes in 
keyaki and sugi wood could not be determined. However, 
the temperature dependency may have been due to the 
extractives responsible for the characteristic color of wood, 
which are usually expressed as chromatic parameters such 
as a* and [29]. Characteristic extractives possibly give 
different E a to different wood species and to different color 
parameters such as L*, a*, /?*, and A E* b . 

From the results of this study, we can conclude that 
color changes in keyaki and sugi wood during heat treat¬ 
ment can be successfully analyzed using a kinetic approach 
combined with the TTSP method. This strategy allows 
prediction of color changes at any of the temperatures and 
durations used in the present study using E a and a T Tem¬ 
perature dependencies expressed by b T must be considered 
when needed. 


Fig. 8 Relationship between 
vertical shift factor b T and 
treatment temperature for <3* 
and b*. a Keyaki wood and 
b sugi wood 
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Conclusions 

Color changes in keyaki and sugi wood during heat treat¬ 
ment at 90-180 °C were investigated in this study. The 
values of L* and A E* b decreased and increased, respec¬ 
tively, in both wood species at all treatment temperatures, 
whereas the values of a* and Z?* showed species- and 
temperature-dependent behaviors. Kinetic analysis was 
successfully applied to describe these color changes. Cal¬ 
culating the values of apparent activation energy and shift 
factors for time-temperature superposition allowed pre¬ 
diction of color changes in the present experimental con¬ 
ditions. The present study revealed and quantitatively 
evaluated differences in color change behavior between 
species and temperatures. The results provide a more 
comprehensive understanding of color changes in wood as 
a function of time and temperature. 
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